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TOMATO YELLOW LEAF CURL VIRUS 
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ABSTRACT. Mathematical model has been developed and analysed for the interaction be- 
tween tomato yellow leaf curl virus (TYLCV) and tomato plants under the influence of 
Bemisia tabaci. Positivity and boundedness of the solution has been checked to ensure our 
model being well posed and then we computed the basic reproductive number Ro using the 
next generation matrix method. Also, both local and global stability analysis at the disease 
free equilibrium points of the model has been done. By constructing suitable Lyapunov 
functional and using LaSalle’s invariance principle, global stability of endemic disease equi- 
librium was obtained. The results show that the disease free equilibrium point (DFE) will 
be both locally and globally asymptotically stable if Ro < 1 but unstable if Ro > 1 and 
the endemic equilibrium point (EE) will be globally asymptotically stable if Ro > 1 and 
unstable if Ro < 1. Finally some numerical simulations were done to validate our theoretical 
outcomes, and the epidemiological implications of the key outcomes were briefly discussed 
in the last section. 


1. INTRODUCTION 


Tomato yellow leaf curl virus (TYLCV) is the name given to complex virus species that 
causes disease known as tomato yellow leaf curl disease (TYLCD). Tomato yellow leaf curl 
virus (TYLCV) is a major threat to the tomato crop in many tropical and sub-tropical 
regions|[12]. TYLCV is transmitted by the insect vector Bemisia tabaci (whitefly). The virus 
accounts for huge losses in quantity and quality of tomatoes if unchecked. Incidences as high 
as 100% with undesirable consequences of crop failure have been recorded[15, 5]. The virus 
was detected first in Israel around 1930, and currently it affects about 30 countries around 
the world that grow tomatoes. Symptoms of TYLCV that exhibited by tomato plants is 
reduction in leaf size, leaf curling upward, severe stunting and distortion associated with 
interveinal chlorosis, observed mainly on the upper portion of plants|16, 13].The disease is 
among the major virus diseases that causes low yields of tomatoes in Tanzania, especially, 
in farmer’s fields. Disease occurrence of about 100% had been reported in some regions in 
Tanzania mainland([9, 14]. 
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2. FORMULATION OF THE TYLCD MODEL 


Plant population has been divided into three categories in SEI type structure, reflecting 
the disease status: susceptible (healthy), exposed (latent infection) and infectious. Similarly 
we have defined categories for the vector population in SI type structure as susceptible 
(non infectious) and infective which links the epidemiology of disease with the population 
dynamics of the vector. In SEI, the total tomato population Nr was subdivided into three 
sub-populations; tomato plants that are susceptible to infection with TYCLV Sr, those 
exposed to TYCLV Er and infectious tomato plants Ir. That is Nr = Sr+Er+/Ir. The total 
whitefly vector population Ny is sub-divided into two sub-populations; susceptible whitefly 
vector population Sy and infectious whitefly vector population Iy. That is Ny = Sy + Iv. 
The transmission dynamics of TYLCV is summarised in the compartmental diagram in 
Figure 2.1. 
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FIGURE 1. Compartmental diagram for the transmission dynamics of TYLCV 


We assume that healthy tomato plants are planted/replanted at the rate A . Plants either 
die naturally at a rate ø or move to the exposed class after acquiring TYLCV through the 
contact with infectious whitefly vector at a rate 6,. The exposed tomato plants are either die 
naturally at rate ø or move to the infectious class at rate G2. Whitefly vectors are recruited 
by two different ways that is by birth and immigration. Whitefly vectors born at rate of A 
and immigration at rate of r. Probability that whitefly vector born susceptible is b while 
probability of being born infected will be 1—b. Probability that whitefly vector immigrating 
susceptible is p while probability of immigrating infected will be 1 — p. They either die 
naturally at a rate y or emigrate at rate e. They move to the infectious class after acquiring 
TYLCV from the infected tomato plants at a rate n. 
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The model assumes that a year-round vegetable production system is used. Replanting took 
place at a rate which exactly balanced those plants removed. This means that we maintain 
a host population at a constant size by balancing the birth and death (or replanting and 
mortality) rates. Farmers plant only healthy varieties of tomato in a garden of carrying 
capacity K, no death of tomato plants before harvesting and the whitefly vectors are assumed 
to remain infectious once they acquire the virus. The variables and the parameters are 
summarised in Tables 2.1 and 2.2. 


TABLE 1. Model Variables 








Variable | Description 

Sr Susceptible tomato plants 
Er Exposed tomato plants 
Ir Infectious plants 

Nr Total plant population 
Sy Susceptible whitefly 

Iy Infectious whitefly 

Ny Total whitefly population 








Table 2: Model Parameters 





Parameter Description 





Planting/replanting rate of tomato 

Total whitefly population birth rate 

Total whitefly population immigration rate 
Probability of susceptible offspring being born 
Probability of susceptible whitefly vector immigrate 
Total whitefly population emigration rate 

Natural death rate of tomato 

Natural death rate of whitefly vector 

Rate of tomato plants to move to exposed class 
Rate of exposed tomato to move to infectious class 
Rate of susceptible whitefly to move to the infectious class 


ZS PPSD m > 


Garden carrying capacity 





From the compartmental diagram in Figure 2.1 and basing on assumptions and descrip- 
tions made we derive the following equations: 
Tomato plants equations 


4 NGALYA AND KUZNETSOV 








dSp Nr, PilvSr 
“ae. = A(1 K ) Nr aSr, 
dE IVS 
(1) FE = EE bab — 0B 
dI 
aE = BoEr = ol rp. 


Vector equations 


I 
woy =A+pt— rey — pSy — ESV, 
(2) dt Ny 
dly nip Sy 
— =(1-—D))A 1— —— — oly = ely. 
Fe = (1 A+ (1 pyr + TEY — ply -ely 


3. BASIC PROPERTIES OF THE MODEL 


3.1. Positivity. In order the TYLCD model to be meaningful epidemiologically, we must 
show that all state variables are non-negative for all times t > 0. We find solution of each 
equation in systems (1) and (2) in their patches for testing positivity. 

Lemma 3.1. Let the initial data be {(Sr(0), Er(0), Ir(0), Sv(0), Iv (0)) > 0} € Q; then 
solution set { S(t), E(t), Ir(t), Sy(t), Iv (t)} of the model system is positive Vt > 0. 
Proof. From the first equation of the tomato plant population of the model system (1), 








dSr Nr, BilyvSr 
He = A(1 K ) Nr aSr. 
dSr Byly 
(3) ore ( Na )Sr(t) 


By separating variables and integrating both sides of the equation (3) above we obtain, 
BI 
S(t) > Sr(0je Tr ~ 7 > 0. Therefore Sr(t) > 0. 


Following similar procedure we can compute all variables Ær(t), Ir(t), Sy(t) and Iy (t) and 
establish that Er(t), Ir(t), Sv (t) and Iy (t) > 0. 
Therefore the solution set {Sr(t), Er(t), Ir(t), Sv (t), Iv(t)} of the model is positive Vt > 0. 


3.2. Invariant region. Since the system is modelling tomato plants and whitefly vector 
populations, we assumed that the state variables and parameters of the model are completely 
non-negative t > 0. The TYLCD transmission model has two compartments where each 
population is treated separately. 

Lemma 3.2. All forward solutions of the of the TYLCD model in R3 enter the invariant 
region Q = Or+Qy where Qr = {(Sr, Er, Ir) €: (Sr + Er + Ir = Nr)}, Qv = {(Sv, Iv) E: 
(Sy + Iv = Ny)}, and Q is the invariant region of the whole system. 

Proof. Let Nr the total population of the Tomato plants Nr = Sp + Er + Ir. Addition 
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of equations in system (1) yield ^r = \(1— *£) — oNr. By separating variables in this 
equation, applying an integrating factor technique we compute the differential equation and 
obtain, 
A+ KX EE A 
Nr = K + C (KTA 
A+ KA A+ KA, (aa, 
when t —> co Nr(t) — “484, fort =0, Nr(t) = Nr(0) which implies that Nr(t) > 0. 


Let Nr be the total population of the whitefly vector, then Ny = Sr + Iy. Addition of 
N 
equations in system (2) yield JF = 2(bA + pr — 1) — (y + £). By separating variables in this 


equation, applying an integrating factor technique we compute the differential equation and 











obtain, 
Ny = 20A tpr- 1) | getet) 
pre 
2(bA + pr —1 ZONE pT — 1), iope 
eet me ) + (ny (0) — % E Yet, where (p+e) £0. 


when t — œ Ny(t) — 2A?) for t = 0, Ny(t) = Ny(0) which implies that 
Ny(t) > 0. 

Therefore since Nr(t) > 0 and Ny(t) > 0 , hence the set {(Sr, Er, Ir) € R3}; ((Sv, Iy) € 
RÊ? )} is a positive invariant set in Q. 


4. MODEL ANALYSIS 

















dI dI 
4.1. Disease free equilibrium point D°. At disease free, di a oe 0, 
dS 
Sa = \(1— 42) — S22 — o Sp, (1 — FE) — PP — o Sp = 0. 
But since at disease free, Iy = 0, Sr = Nr therefore A(1 — Sr) — es = 0. Making Sr the 
subject, we obtain Sp = aK +x: Considering following equations, a = A4 ox urey = 
dI 
Sy — ESy, a= (1—bd)A+(1—p)rt+ wey — ply — ely. 
Irs 
(4) Akp IE euis 
Ny 
Irs 
(5) (1-0)A+ (1 -pr + TEY — ply — ely =0. 
V 


Addition of these equations (4) and (5), results bA + pr — pSy —eSy + (1 —b)A + (1 —p)r — 
ply —ely = 0. But Iy = 0. 
bA +pT-— pSy —ESy +(1—b)A+(1—p)r =0, Sy = HE, 


Therefore the D? of the model system is given by (A 0,0, ee 0). 
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4.2. Endemic equilibrium point D*. In the presence of TYLCD, Er 4 0, Ir 4 0, Iy # 
0 our model has an equilibrium point called endemic equilibrium point denoted D* = 
(Si, Ez, 17, 5}, To) 4 0. D* is the steady state solution where TYLCD persist in the pop- 
ulation of tomato plants. For the existence of D*, the elements must satisfy; S% > 0, £7 > 
0, I} > 0, S% > 0, > 0. We find the endemic equilibrium point by setting the right side of 


the model system equations (1) and (2) equal to zero. Thus; 
[I= Bao(bA+pr)— (pte) 
Bam- EL) (820) 


* 





, 








* ol*T 
Ep Bo ? 
x _  (bA+tpr) Nọ 


T nlp t(pteyNy? 
qe = NEUE+(e+e))(G-b)A+(1-p)r) +ITNV (b+) 
Br (e+e) (nlp+(ete) Ni) : : 
For a positive endemic equilibrium point, the conditions 4(1 — Nr) ao (B2 + 0), Boo (bA + 


K Be 
pt) > (p +e) and BgnX(1 — NF) > (B2 + 0o) must hold. 











4.3. The basic reproduction number Ro. The basic reproduction number (Ro) is the 
expected number of secondary cases produced by a single infection into a completely sus- 
ceptible population[1]. [3] defined Ro as a spectral radius (o(FV7!)) of the next generation 
matrix. This is a parameter which usually tells if the epidemiologic might persist or die out 
within a purely susceptible population. The epidemiological criterion of Ro is that if Ro > 1 
then the disease free equilibrium point is unstable and can invade the population and persist 
for long time [17]. There are several methods for calculating Ro but in this study the next 
generation operator method as proposed by [17] is used. 

The disease transmission model comprised of the system of equations X} = Fi(x) — V;(x), 
where F;(x) is for new infection and V;(a) is for remaining transfer terms. To obtain the 
matrix F and V (Jacobian matrices), we then find derivatives of vectors Fj(x) and V;(x) 
respectively. 

From the model equation, 


Bilv ST 
Nr 


F;(x0) = 0 
(1—D)A+(1—p)r + EY 


The Jacobian matrix at D° is 
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From the model equation, 








boEr + o Er 
V}(2)= | oly — b2Er 
ply + ely 
The Jacobian matrix at D° is 
b2 +0o 0 0 AT 0 0 
V=| -b6 o 0 |, andV7= | 2 0 
1 


Multiplying F and V~! yields a next generation matrix FV~! as shown below, 


B 
0 0 ar 
Fv '=| 0 0 0 
Ban n 
DD o 0 


Computing the maximum eigen value of the next generation matrix we get Ro = ,/ a 
Ro is a threshold parameter that indicates the average number of infected vectors and in- 
fected hosts caused by a cross-infection of one tomato plant host or one whitefly vector when 
the other population consist of only susceptible [10]. The square root arises from the fact 
that dual generations are necessary for transmission of TYLCD to take place, i.e. From an 
infectious tomato plant to a susceptible whitefly vector and then from an infectious whitefly 
vector to susceptible tomato plant(host) [10]. 


4.4. Stability analysis. 


4.4.1. Local stability of the disease free equilibrium point. We examine the stability of the 
AK bA 
K+X? 0, 0, p 

o p 








disease free equilibrium point D° = ( 
Kinene et al. (2015). 
Lemma 4.1. The disease free equilibrium D? is locally asymptotically stable if Ro < 1 but 





br, 0) by employing the method used by 


unstable if Ro > 1, where Ro is the basic reproduction number. 

Proof. We linearize the model system (1) and (2) by computing its Jacobian matrix Jpo. 
The Jacobian matrix is computed at disease free equilibrium point by differentiating each 
equation in the system with respect to the state variables Sr, Er, Ir, Sy, and Iy. Therefore 


the Jacobian matrix is 


eee 0 0 -bı 
0 —b2—-0 0 0 pı 
(6) Jpo = 0 By =o 0 0 
0 0 Si Sa re TE] 0 
0 0 n 0 —p-—€é 
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It is clear from equation (6) that the first and second eigen values are Ay = —o and 
Ao = — (ee) + +€). This can be reduced to a (3 x 3) matrix 
-2—0 0 bı 
(7) Jpo = Bo —0 0 
0 n -pr-e 


Its characteristics polynomial is 
(E+ p+ 20+ Bo)A? + (0? + (2u + 2€ + Bo)o + Bole + m))A+ (E+ u)o? + (e + 4) B20 — Bi Bon + 
(e+ wt 2a + Bo)d? + (0? + (2u + 2e+ Bolo + Bole + w))A+ (E+ po? + (e + u) Boo — 81 bon. 
a, = E+u+20+b2, ay = (07 +(2u+2e+f2)o+fo(etp)), a3 = (e+u)o?’+(e+u)b20— b1 Ban. 
Consider (€ + u + 20 + 62) > 0. Since all model parameters are positive, it is obvious that 





























a, is positive (a, > 0). 

Consider (o? + (2u + 2€ + B2)o + Go(e + u)) > 0. From the point that all model parameters 
are positive, it is obvious that az is positive (az > 0). 

In order that az being positive, (£ + u)o? + (e + 1) B20 — B1Bon > 0. 








(e + [jo + (e at [t) Boo > 6 Bon, 1> EE 1> m 
a1az—as = (e+u +20 +o) * (0? +(2u+2e4+ B2)o+ Bo(et+p)) —((e+p)o? + (E+) G20 — 81890). 
After simplification yields a1a2 — a3 > 0. 





By Routh-Hurwitz criteria all eigenvalues have negative real parts if Ro < 1 thus making 
the disease free equilibrium locally asymptotically stable. 


4.4.2. Global stability of the disease free equilibrium point. We examine global stability of 
the disease free equilibrium using the theorem proposed by [2] similar as in [10]. Thus we 
rewrite our model and mention two conditions, if satisfied will guarantee global asymptotic 
stability of the disease free equilibrium. 


dX 
“y eo F(X, Z), 

(8) F(X,Z)= 4 # 
GF =X), G(X,Z)=0, 


where X = (Sp, Sy) € R? denote susceptible populations and Z = (Er, Ir, Iy) € R? 
denote infected populations. D° = (X*,0) represent the disease free equilibrium of this 
system. The conditions (i) and (ii) below guarantee global asymptotic stability: 

dX 
(i): for ao F(X,0), X* is global asymptotically stable. 

(ii): G(X, Z) = DzG(X*, Z)Z — G(X,Z), Ĝ(X,Z) > 0 for, (X,Z) €Q, 
where DzG(X*, Z) is Metzler Matrix (the off diagonal elements are non negative) and also 
the Jacobian of G(X, Z) taken in (Er, Ir, Iv) and evaluated at (X*,0) = (A HE, 0,0,0). 
If the system satisfies the above conditions, then according to [2] the following theorem holds. 
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Lemma 4.2. The equilibrium point D? = (X*,0) of the system (8) is globally asymptoti- 
cally stable if Ro < 1 and the conditions (i) and (ii) are satisfied. 
Proof. We now start our proof by describing new variables and breaking the system into 
subsystems. X = (Sr, Sy) and Z = (Er, Ir, Ivy). From equation (8) we have two vector 
valued functions G(X.Z) and F(X,Z) given by: 


AG- Oe = ayer =e 
bA + pr — ey — Sv — Sy, 


alse- Er — Br, 
G(X, Z) = BaBy —olr, 
(1-b) A+(1-—p)r + ue — ply — ely. 


7 
For condition (i) we considering the reduced system —— = F'(X,0). 


dt 
dS N Ivs 
—— =X(1 — 4) — 88" — Sr, 
(9) F(X,Z)= 4 ff, ae 
Hw At er Tt — Psy — ESy. 








xX* = ( AK bA+pr 


RIN eee ) is globally asymptotically stable equilibrium. To prove this, we find 


dX 
the solution of ape F(X,0). When we find solution of the first equation from (9) yields 


Sr = ea + (Sr(0) = xe Oo 
bA 
the second equation of (9) yields Sy = “er ae S(O) = CARED eer 


pre 
bA+pt 
=a ift — oo regardless of the initial condition. Thus, is globally asymptotically stable. 


KEA 








t that approaches Ak oy ft — oœ. panes solving 








t that approach- 





es 


Meanwhile for condition (ii) we compute G(X, Z) = DzG(X*, Z)Z — G(X, Z), and then 
prove that Ĝ(X, Z) > 0, for (X,Z)€Q. 


Anse — PoE, — o Er 
G(X, Z) = "Brali 
(1-—b)A+(1-—p)r+ wey — ply — ely 


At disease free equilibrium point, 


—(ßı +o) 0 By (01+ 0)Er + Bly 
DzG(X*,0)= Bo —O 0 5 DzG(X*,0)Z= BoE — olr 
0 n —ly+e) nIr— (y + £)Iv 


By applying the formula G(X, Z) = DzG(X*, Z)Z — G(X, Z), 
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-(bı + o)Er + Bly AT — pEr — o Er 
G(X, Z) T BoE E olr Ez] ae = olr 
nlp — (ẹ + £€)Iv (1-bA+(1-p)r + 2 — gly - ely 
; ily (1 — 22) 
CA A= 0 
m: ((1—6)A + (1—p)r + MSY) 
Since (1— xr) > 0 then 6,Jy(1- =) > 0, meanwhile for nlp > ((1—b)A+(1—p)r+ wey) 


then nlp — G- b)A+(1—p)r+ lps) > 0. Therefore G(X, Z) > 0 and hence this completes 
the proof of lemma 4.2. 


4.4.3. Global stability of endemic equilibrium point. Global stability of the endemic equi- 
librium point D* is analyzed by constructing a suitable Lyapunov function. We prove 
for global stability of endemic equilibrium point following the approach used by [8] and 
other several epidemiological models. We consider the Lyapunov function of the form 
L => Pi(Ri — RF In(R;)). 

So we define the Lyapunov function as L(Sr, Ir, Sv, Iv) = Pi(Sr — Sž ln(Sr)) + Po(Er — 
Ex \n(Er)) + Ps(Ir — I% ln(Ir)) + Pa(Sy — S$ In(Sy)) + Psy — Tf InUy)). 


We can now verify the condition a <0. 
dL dS dE dI 
Proof. ie derivative ce is ae P\(1- aay +P(1— zi) =a +P(1— ers + 
ra pare ee Sew E ee 
a-pa BS 


= P- ZAU- 4) - ox — oSr) + Pa(1 — F(E — pEr — oEr) + P3(1 — 
1) (BB — ofp) + Pi(1 — St)(A + pr — UES” — ySy — eSy) + Pal EEA 
p)t + ue — ply — ely). 

At an endemic equilibrium point D* we have \(1— 42) = (Ax +0) Sh, Vr = (Ba +0) Et, 











K 
bB2Er. = oI}, bA + pr = (ev +y+e)S% and ((1—b)A+(1—-p)r4 ursy = = (y + e). 
= P(1- 3)((A¥ +0)S} — (AM + 0)Sr) + Pa(1 — F) ((62 + 0) Bf — (b2 +0) Br) + P;(1 — 
E)(olp—olr)+Pa(1— ge) (HE ae (HE +p-+e)Sy+Ps(1—72)((pte)ly—(pte)lv). 





= —P,(1— £)?o)Sp— P(e (1 = se) she) — P2(1— $£)?(Bp +0) Er — Ps(1— £)?o Tp — 
Pi(1 — 3¢)°(p + )Sy — preseg — =)(1— 4) P1- Evle +e). 

= —P,(1— 4)°o) Sp — Pall HG +0) Ep — Pa(l — Poly — Pill- Bp + )Sy - 
P3(1— Ty Ee) P (Sx, Bas Pe, Sva V, 

where P(Sr, Br, Ir, Sv, Iv) = (—P, (4 (1— 3£)(1— 34 *) — PME (1 4)(1- 4) < 0. 
The function P(Sr, Er, Ir, Sv, Iv) is non-positive by considering the method used by [8]. 
Thus P(Sr, Er, Ir, Sv, Iv) < 0 for all Sr, Er, Ir, Sy, Iy > 0. Hence, < 0 and = = 
if and only if Sp = Sp, Er = Eñ, Ir = Iž, Syv = Sy and Iy = Jy. Therefore the largest 
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invariant set such that ae = 0 is the singleton {D*} which is our endemic equilibrium point. 
By LaSalle’s invariant principle in [11] we conclude that D* is globally asymptotically stable 
if Ro > 1 and unstable if Rọ < 1. So, we establish the following lemma. 

Lemma 4.3. The endemic equilibrium point D* of the model systems (1) and (2) is globally 
asymptotically stable if Ro > 1 and unstable otherwise. 


5. NUMERICAL SIMULATION AND DISCUSSION 


The main intention of this paper is to assess the impact of number of bemisia tabaci 
interacting with tomato plants in transmission of TYLV. In order to support the analytical 
results, numerical simulations were performed using different values for the initial number of 
vectors but constant number of host population (tomato plants) and using numerical value 
of parameter taken from literature while some were estimated. Fixed variables used are 
Sr=2000, Er=I7=0, while other variable were varying (Sy=Iy=500, 2500 and 5000). The 
final time is t;=—90 days, and the rest parameters values are shown in the Table 5.1. 


Table 3: Parameter values for the model 














Parameter Value day! Reference 
À 0.01 4 
0.0118 18] 
E 0.3 4 
b 0.5 T 
p 0.8 4 
T 0.267 4 
p 0.01 y 
o 0.015 10] 
bi 0.01 4 
Bo 0.075 
n 0.33 Estimated 
K 2000 Estimated 





Graphs of susceptible tomatoes, exposed tomatoes and infected tomatoes were plotted a- 
gainst time by varying number of vector in the field. Computations were done using Matlab 
with the ode45. This function implements a Runge-Kutta method with a variable time step. 
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FIGURE 2. Graph of susceptible tomato plants vs. time for the first 90 days of growth 
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FIGURE 3. Graph of exposed tomato plants vs. time for the first 90 days of growth 
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INFECTED TOMATO PLANTS 
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FIGURE 4. Graph of infected tomato plants vs. time for the first 90 days of growth 


Numerical simulation results in Figure 5.1 shows significant decrease in the number of sus- 
ceptible plants due to increase of number of vectors in the field. Figure 5.2 shows an increase 
in the number of exposed plants due to increase of number of vectors in the field and Figure 
5.3 shows also the increase in the number of infected plants due to increase of number of 
vectors. This implies that the increase of number of vectors around the field lead to decrease 
of amount of tomato yield in that field due to TYLCV. 


6. CONCLUSION 


We computed for the basic reproductive number (Rọ) that can be used to determine 
whether eradication of TYLCV will occur (Rp < 1) or if disease will continue to increase 
(Ro > 1). Numerical simulation of this model allows us to describe numerically depen- 
dence between the number of vectors and infection transmission to tomato plants. A better 
understanding of plant-virus-vector systems was obtained with the development of model 
incorporating crop fields and vector populations with migration of vectors between fields. 
We notice that in order to lessen the infections, intervention strategies need to be focused 
on killing vectors so that we reduce the contact between bemisia tabaci and plant host. 
Disease controlling in one tomato field may impact the population dynamics and propor- 
tions of susceptible and infectious vectors and affects disease dynamics in other fields due to 
vector migrations |6]. There is a need of effective insecticides application that would reduce 
the whitefly population and keep the tomato population stable. Therefore, this modelling 
approach can be modified and then used to evaluate control strategies of the disease. 
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